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Abstract 

In the limit of a high temperature T and a large quark-mass M, implying a small gauge coupling 
g, the heavy quark contribution to the spectral function of the electromagnetic current can be 
computed systematically in the weak-coupling expansion. We argue that the scale hierarchy 
relevant for addressing the disappearance ("melting") of the resonance peak from the spectral 
function reads M S> T > g 2 M > gT ^> g 4 M, and review how the heavy scales can be integrated 
out one-by-one, to construct a set of effective field theories describing the low-energy dynamics. 
The parametric behaviour of the melting temperature in the weak-coupling limit is specified. 
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1. Introduction 

Given possible applications in heavy ion collision experiments [1] , the last 20 years have 
seen a huge amount of phenomenological work on the properties of heavy quarkonium at 
temperatures just above the deconfinement crossover (for a recent review, see ref. [2]). 
During the last couple of years, these phenomenological works have been complemented 
by more theoretical investigations, aiming at a QCD-based approach to the problem. The 
purpose of this talk is to distill a basic message from some of the latter studies |3J-|9 r 

Let us start by defining the observables we are interested in. The heavy quark contri- 
bution to the spectral function of the electromagnetic current can be written as 



oo 

PV {Q) = JdtJ d 3 x e «** /i[XaO,i M (0)]\ 



where = tp 7 M -0; is the heavy quark field operator in the Heisenberg picture; (...) = 
Z^ 1 Tr[(...)e~^ H ] is the thermal expectation value; f3 = 1/T is the inverse temperature; 
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Fig. 1. Phenomenological results for (left) the spectral function py from ref. [4]; and (right) the dilcpton 
production rate dN^- M + /d 4 xd 4 Q from ref. [6], In both cases Q = (ui, 0). 

and we assume the metric convention (H ). This spectral function determines the 

production rate of muon-antimuon pairs from the system |10j . 
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where Z is the heavy quark electric charge in units of e, and tib is the Bose-Einstein 
distribution function. In the following we assume, largely for notational simplicity, that 
the muon-antimuon pair is at rest with respect to the thermal medium, i.e. Q = (u>, 0)P1 
Furthermore, M denotes the heavy quark (charm, bottom) pole mass. 

The ultimate goal of the study would then be to compute pv and dJV„-„+/d xd Q 
in a certain energy range around the two-particle threshold, say uj ~ (1 . . . 3)M, and in 
a certain temperature range, say T ~ (0.1 ... 1.0) GeV. The phenomenological interest 
of the problem lies in the fact that the observables mentioned are believed to undergo a 
qualitative change in this temperature range, thus constituting a sensitive "thermometer" 
for the system; examples of possible patterns are shown in Fig. [TJ 

Now, given that we are interested in fairly low temperatures and that the charm quark 
mass is fairly moderate compared with the QCD scale, it would be interesting to tackle 
the problem with lattice techniques (see, e.g., refs. [H]). Unfortunately, lattice methods 
are Euclidean, and even for a perfectly known Euclidean correlator Gy(r), r S [0, 0], it 
is strictly speaking not possible to invert the relation to py, viz. 

~'dcu , cosh(f-r)c 
- PV{U) sinh^ 



OC 

Gv(r) = f - 



without inserting further input (for instance, in momentum space, one must set e I/3w " = 1 
in order to get the correct analytic continuation from uj n to a;). Another possible strong 



1 For a non-zero total spatial momentum q, with < |q| <C M, the main modification would be a shift 
of the two-particle threshold, from ui fs 2M to U) fa 2M + q 2 /4Af (see, e.g., ref. [6]). 



2 



coupling method, an analysis through a gravity dual (see, e.g., ref. [12j). is not available 
for QCD. In the following we therefore resort to the weak-coupling expansion, which at 
least is a theoretically consistent framework; in fact, particularly for the bottomonium 
system, experience from zero temperature (see, e.g., ref. [13]) suggests that numerical 
convergence could be reasonable as well. 



2. Scale hierarchy 



A systematic weak-coupling computation can be carried out once the scales that en- 
ter the problem have been identified. In the heavy quark-mass and high-temperature 
limit, the QCD coupling constant g is small, and the quarkonium states resemble non- 
relativistic systems such a positronium. Thereby we treat the heavy quark mass M as the 
"hard" scale; the inverse Bohr radius, or relative momentum Mv ~ g 2 M as the "soft" 
scale; and the binding energy Mv 2 ~ g A M as the "ultrasoft" scale. We assume the QCD 
scale to be at most as large as the ultrasoft scale. 

The energy scale with which we probe the system, u>, could in principle be anything. 
However, conceivably the most interesting range is \u> — 2M\ ~ g 4 M, i.e. deviations from 
the two-particle threshold by at most the binding energy; this is where the quarkonium 
resonance is to be found. In the following, we concentrate on this regime. 

Finally, we need to decide how temperature relates to the zero-temperature scales. 
Maybe quarkonium melts when T ~ g 2 M, given that then hard particles from the plasma 
have enough momentum to kick constituents out of the bound state? However, such a 
temperature cannot be high enough, since kicking happens through interactions, and 
those are suppressed by g. On the other hand, increasing the temperature just a little 
bit, to gT ~ g 2 M, is certainly enough to dissociate the bound state, given that then even 
the Coulomb binding gets Debye screened. In fact, removing completely the Coulomb 
binding is an overkill, since a finite width could arise already earlier on. To summarise, 
the relevant range is somewhere between these two limits, g 2 M < T < gM [6]. Therefore 
we now assume the equivalent hierarchy 

M > T > g 2 M > gT > g 4 M , 

and integrate out the hard scales one-by-one. 



3. Integrating out the scale M 



Integrating out the hard scale M yields an effective theory called NRQCD [14 15]. The 
topic is very well developed by now, and we only list some basic steps here; many details 
and references can be found, e.g., in the review [16) . 

Starting from the QCD Lagrangian, £qcd = "^{i^D^ — M)ip + £ g i U on, one can first 
carry out a Foldy-Wouthuysen transformation: 




where rlj = dj+igAj, Dj = dj—igAj, and 9,<p are two-component spinors. Expanding 
in 1/M and dropping total derivatives, this yields 
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The heavy masses on the first row can be shifted away by — » e~ lMt 6, <f> — > e lM '</>. 
The mixed terms on the second row become rapidly oscillating after this shift, e 2tMt 9^a ■ 
gE (f> — e~ 2tMt <ft<j ■ gE 9 , and describe hard reactions; such oscillations can be integrated 
out yielding four-fcrmion operators. Finally, at the loop level, we need to correct for 
the local expansion through matching coefficients, which concerns both parameters, e.g. 
M = m^g(m^g) (l + i g 2 C , F/47r 2 + ...), as well as composite operators, e.g. 



4. Integrating out the scale T 
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The next step is to integrate out the hard thermal scale T. In the heavy-quark sector, 
this does nothing dramatic, just corrects the parameters. In fact, since the thermal loops 
can be computed within NRQCD, all effects are power-suppressed in mass, e.g. jTT] 

12M 

In the gauge field sector, the same step yields nothing but the Hard Thermal Loop 
effective theory [Tgj: £ g i UO n — ► £gluon + £htl + •••• 

5. Integrating out the scale g 2 M 

The crux of the problem is the integration out of the scale g 2 M, yielding an effective 
theory called pNRQCD [19 20 . (To be more precise, all gluons of energy or momentum 
~ g 2 M are integrated out, while the composite quarkonium fields left over as dynamical 
degrees of freedom can still have relative momentum ~ g 2 M; it is only their energy that 
is small, ~ g 4 M.) Although well established, this step is in many ways more delicate than 
the previous ones; we only outline the basic ideas here (for reviews, see refs. [21116) ). 

The pNRQCD setup is special in that in only applies to particular Green's functions; 
in the thermal context, we can say that the thermal expectation value is restricted to the 
quark-antiquark sector of the Fock space, Tr[(...) exp(— [3H)] — > Tr[P 3 » 8) 3(...) exp(— f3H)]. 
Thereby exponentially small effects ~ exp(— (3M) are omitted. For the quark-antiquark 
pair, were represent colour as 3* ® 3 = 1 © 8, and space coordinates as Xi = X + r/2, 
X2 = X — r/2. The singlet field, S, is diagonal in colour-space, S = S 1n c xn c / V^c, while 
the octet field is traceless, O = O a T a /\/Tf, where Tp defines the generator normalization 
through Tr[T a T h ] = TpS ab . The interpolating operators for the component fields S,O a 
can be chosen as 

t (x 2 ,i)W X2iXl 0(x 1 ,t) ~ Z a *(r)S(X,r,t) , 
e t (x 2 ,t)W/ X2 , x T a IFx,x 1 0(xi,t) ~ zl(r)O a (X,r,t) , 
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where W is a straight Wilson line. Subsequently a general Lagrangian is written for the 
fields S, O as well as the ultrasoft gauge fields, respecting gauge invariance and expanded 
in 1/M, l/g 2 M ~ r: 



C ~ Tr I S f 



V 2 V s {1 \r) 



M M 



V A (r) [Qtr-sES + Str-gEO] 



S + Qt 



V B (r) 



M M 



[O f r -.gEO + O^r • .gE] \ + 



O 



1 



where T>q is the covariant derivative in the adjoint representation, i.e. iT>qO = idoO + 
g[Ao(X,t), O]. As indicated, gauge fields only depend on the center-of-mass coordinates 
X, t, whereas dependence on r can be Taylor-expanded. As far as the structure on the 
second row is concerned, we remark that at leading order Va — Vb = 1, as can be 
inferred from a tree-level matching of a three-object vertex to NRQCD, on which side 
the covariant derivative can be expanded as iDo — > iDo + r • gE + ... . 

Now, the effective theory contains many "potentials", V s , Vo , Va, V Bl which are 
to be treated as matching coefficients. They can be computed from various Wilson loops 
with t ^> r (since dt ~ g A M <C |V| ~ g 2 M), i.e. in the static limit. We return to the 
determination of Vs presently. 

For the computation of the spectral function py , a representation of the vector current 
within pNRQCD is also needed. The vector current is a local object, corresponding to 
r = 0. At zero distance, Z s (0) = N c , and the relation between pNRQCD and NRQCD 
fields becomes unambiguous. In addition, spin indices, which have been suppressed in 
the discussion above, can be added in a trivial way. Thereby the spectral function can 
indeed be computed within pNRQCD (for a concise yet explicit presentation at zero 
temperature, see ref. [2"2"]). 



6. Integrating out the scale gT 

Even though our ultimate goal is to treat the situation g 2 M > gT, in which case the 
scale gT is integrated out within pNRQCD, it will be illuminating to start by considering 
the case g 2 M ~ gT. Indeed the scale gT can then be integrated out together with the 
scale g 2 M , and the problem boils down to determining the potentials of the previous 
section in the presence of non-zero Debye screening (since m D r ~ 1). 

Luckily, this problem was considered already in ref. [3]. In that paper the result was 
called a "real-time static potential" , and was in general denoted by V> (t, r); the potential 
Vs°\r) equals the infinite-time limit thereof, vj°\r) = lim^oo V>(i, r). Thereby [5] 
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where to d ~ gT is the Debye mass, Cf = 4/3, and 
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Let us comment, in passing, on the existence of an imaginary part in the potential. 
For r — > oo, the real part becomes a constant, corresponding to twice the thermal heavy 
quark mass shift [23] ; in complete analogy, the imaginary part at r — > oo equals [7] twice 
the known heavy quark width [24j . In some sense, the physics of the imaginary part is 
also closely related to (single) heavy quark energy loss / drag force / diffusion coefficient, 
on which a lot of work has been carried out recently (see, e.g., refs. [55]): in fact the 
essential graph can be depicted in the same way, 



X 



but the integral over the spatial momentum of the gluon is weighed differently. 

What we really wanted to do in this section, however, was to integrate out the scale 
gT for g 2 M > gT . One possibility would now be just to expand the previous potential 
in m D r ~ gT/g 2 M < 1. It is illuminating, however, to really derive the result in two 
separate steps, like the effective theory setup requires; the imaginary part turns out to 
have an 1/e divergence, and it is a nice crosscheck to see that this cancels in the end [9]. 

The first step consists of the computation of the potential Vs but with to d < |p| < T, 
where |p| ~ 1/r ~ g 2 M is the gluon momentum. Thereby we are advised to expand the 
static limit of the time-ordered HTL gluon propagator, 



tf>£>(o,p) 



i 



p 2 +m^ |p|(p 2 + ml,) 2 
to leading non-trivial order in m^/p 2 , and the imaginary part becomes 



iirg Cpm^T 
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where p, 2 = 47r/x 2 e _7E is the scale parameter of the MS scheme. 

In the second step, the scale gT is integrated out within the pNRQCD effective theory. 
The effect comes from a loop where the singlet splits into an octet and a colour-electric 
field. This is mediated by the vertex containing the potential Va (which, as mentioned, 
equals unity at leading order). For the loop momenta relevant for the integration, the 
octet propagator is ~ S(p°), so again the static limit of the (colour-electric) gluon prop- 
agator is needed. This leads to 



6 gT VP(r) 



ing 2 CFm 2 )T 



ing CFm D T 



d 3- 



2c, 



Cp- 



(2irf-^ |p|(p 2 



24tt 2 



In 



ml) 2 
- 2 In 2 



Summing the two parts together, we arrive at [S] 
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Re K C°>(r) = -^i + ..., 
Im^) = -^^M 



4. 3 \ -,„„■ ' 3, 1 

7. Physics at the scale g A M 

Having integrated out all but the ultrasoft scale, let us now inspect the singlet propaga- 
tor at the lowest energies. We can expect a resonance to melt once the width appearing in 
the propagator becomes as large as the binding energy. Of course, this condition cannot 
be posed in an exact way, and for real- world applications we would indeed be interested 
in computing py in some macroscopic energy and temperature range; nevertheless, for 
the purposes of this talk, we restrict to the task of estimating the melting temperature. 

Parametrically (omitting Cf/^ and numerical factors), the equality of the binding 
energy and width from the end of the previous section reads 

a 2 1 

if Irr, 2 2 1 1 

— ~ g 1 m D r in ■ 



Inserting r ~ l/g 2 M and m D ~ gT, we get 
ff 4 M 3 - T 3 In %— . 

For small enough g this can be solved in a leading-logarithmic approximation, yielding 
T~ gi(\n±)-^M . 

This result was first obtained (without the logarithm) in ref. [5] , and it indeed lies within 
the range g 2 M < T < gM discussed in Sec. as must be the case [5J. 

8. Conclusions 



The main purpose of this talk has been to underline the question posed by the title. 
Indeed, it appears that quarkonium physics at high temperatures is important enough 
phenomcnologically to deserve some theoretically-minded consideration as well. 

As far as an answer to the question is concerned, we have tried to obtain one within the 
framework of the weak-coupling expansion. The most important point to realise is that 
a hierarchy can be found between the different dynamical scales affecting the problem. 
In fact, supplementing the discussion so far with the non-perturbative colour-magnetic 
scale g 2 T, the hierarchy relevant for the melting of quarkonium can be written as 

M > T > g 2 M > gT > g 2 T > g 4 M . 

Furthermore, for energies around the two-particle threshold, \u> — 2M\ ~ g 4 M, we 
have argued that the way to exploit the scale hierarchy is to generalize the effective 
theories known as NRQCD and pNRQCD to finite temperatures [819] . As a result of 
such an analysis, it can be argued that quarkonium melts at T ~ g 4 / 3 M [8]. We find it 
comforting that the phenomenon of melting can thus be confirmed model-independently. 
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For practical applications of the effective theory setup, the issue of the convergence of 
the weak-coupling series must be addressed. For this purpose, it would be helpful to use 
the framework to compute higher order corrections to the quarkonium spectral function. 
Also, the influence of the colour-magnetic scale g 2 T needs to be estimated [5]. It appears 
reasonable to expect that as a result of such work, at least the study of the bottomonium 
resonance can ultimately be promoted to a quantitative level. 

Finally, irrespective of the numerical convergence of the weak-coupling expansion, it is 
perhaps useful to emphasize the new qualitative features that this approach has unveiled. 
In particular, the fact that the real-time static potential has an unsuppressed r-dependent 
imaginary part at finite temperatures should in our opinion be taken into account also 
in phenomenological potential model studies. 
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